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ABSTRACT 

Context. Matter leaving the donor during mass transfer spins up the gainer and creates a hot spot in the impact area. If the kinetic 
energy of the enhanced rotation combined with the radiative energy of the hot spot exceeds the binding energy of the system, matter 
can escape from the binary. 

Aims. We calculate the amount of mass lost during eras of fast mass transfer We simulate the distribution of mass-ratios and orbital 
periods for interacting binaries with a B-type primary at birth where mass transfer starts during hydrogen core burning of the donor 
/Wefhods. We used the initial distributions of primary mass, mass-ratio and orbital period established in a previous paper The amount 
of time the binary shows Algol characteristics within different values of mass-ratio and orbital period has been fixed from conservative 
and liberal evolutionary calculations. We use these data to simulate the distribution of mass-ratios and orbital periods of Algols with 
the conservative as well as the liberal model. 

Results. Rapid rotation and hot spots are frequently observed at the surface of the gainer in a semi-detached binary. The mass transfer 
rate for low-mass binaries is never sufficiently large to achieve mass loss from the system. Intermediate-mass binaries blow away a 
large fraction of the transferred mass during short eras of rapid mass transfer 

Conclusions. We compare mass-ratios and orbital periods of Algols obtained by conservative evolution with those obtained by our 
liberal model. We calculate the amount of matter lost according to our model by binaries with an early B-type primary at birth. Since 
binaries with a late B-type primary evolve almost conservatively, the overall distribution of mass-ratios will only yield a few Algols 
more with high mass-ratios than conservative calculations do. Whereas the simulated distribution of orbital periods of Algols fits the 
observations well, the simulated distribution of mass-ratios produces always too few systems with large values. 
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1. Introduction 

Eggleton ( I2OOOI 1 introduced the denomination liberal to dis- 
tinguish binary evolution with mass and subsequent angular 
momentum loss from the conservative case, where no mass 
leaves the system. A catalogue containing nowadays 240 con- 
servative evolutionary tracks can be found at the VUB-website 
(12003-20091) . Van Rensbergen et al. (2008 ) developed a liberal 
scenario in which mass can be lost from a binary during a short 
era of rapid mass transfer soon after the onset of Roche Lobe 
Overflow (RLOF). The joint online catalogue with currently 356 
liberal evolutionary tracks is available at the Centre de Donnees 
Stellaires (CDS). The grid of calculations covers only binaries 
with a B-type primary at birth and initial orbital periods so that 
RLOF starts during hydrogen core burning of the donor: the case 
A of RLOF. Systems undergoing a second era of RLOF after the 
onset of hydrogen burning in the shell of the donor are desig- 
nated as cases A/B. 



2. Initial conditions for binaries with a B-type 
primary 

Van Rensbergen et al. ('2006') used non-evolved systems in 
the 9th Catalogue of Spectroscopic Binaries of Pourbaix et al. 
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(120041 ) to establish the initial conditions for the evolution of bi- 
naries with a B-type primary at birth. We distinguish between 
late B-type primaries in the mass-range [2.5-7] Mq and early 
B-type primaries in the range [7-16.7] Mq. The subscript d is 
used throughout the paper for the binary component which is 
the donor during RLOF, whereas the subscript g is used for its 
mass-gaining companion. The following initial conditions were 
found: 

• future donors follow a normalized IMF: ^(Mrf)=CxM^" with 
Qf=2.35 as given by Salpeter ( I1955I) . 

• initial orbital periods obey a normalized distribution: n(P)=^ 
as given by Popova et al. ( I1982I ). Different distributions 
between small and large initial periods (Van Rensbergen et al. 
2006) do not have to be taken into account since we consider 
only small periods leading to RLOF A in this paper. 

M 

• initial mass-ratios (q=^) follow a normalized distribution: 

^(q)=Cx(l + qY^ ; (5=0.65 for non-evolved binaries with a late 
B primary and 5=1.65 for non-evolved binaries with an early B 
primary, as given by Van Rensbergen et al. ( 120061 ). 

3. Amount of mass driven out of a close binary 

Using a Monte-Carlo simulation starting from a large number of 
non-evolved binaries, the liberal binary evolutionary scenario of 
Van Rensbergen et al. ( 20081 ) has been applied to binaries with 
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a B-type primary at birth and undergoing RLOF A, i.e. during 
hydrogen core burning of the donor The details used in the cal- 
culations can be found in that paper. 

3.1. The contact phase 

A hot spot is created on the gainer's equator or on the edge of 
its accretion disk when the binary is semi-detached. For initial 
mass-ratios around 0.4 or lower (and in some cases for mass- 
ratios as high as 0.6), systems undergoing RLOF A evolve into 
contact soon after the onset of mass transfer, as a result of a rapid 
increase of the outer surface layers of the gainer. The effect of an 
energy stream during contact has been studied by Nakamura & 
Nakamura ( |1987albl l and Packet (119881 available upon request). 
Packet ( 1988) found that the effect has negligible influence on 
the evolution of the binary if the energy stream goes from gainer 
to loser (i.e. if the gainer is the hotter star, which is the case in our 
systems). Hence, we assume that mass transfer continues during 
this phase, without any extra energy transfer 

No hot spot can be created or spin-up generated during con- 
tact since the geometrical impact parameter of the system van- 
ishes in that case. The binary suffers then no mass-loss and does 
not meet the Algol requirements. The contact phase does not last 
long and the system returns into the semi-detached state soon af- 
ter mass-ratio reversal, as a combined result of a decrease of the 
mass transfer rate and the increase of the distance between the 
components. Initial mass-ratios smaller than 0.25 have been ex- 
cluded from our statistical analysis because we assume that con- 
tact in those cases will provoke merging of the system. Initial 
mass-ratios larger than 0.25 undergoing contact have been in- 
cluded for the reasons given in this subsection and also because 
their exclusion would lead to the production of a too small num- 
ber of Algol type binaries. 

3.2. The hot spot 

The accuracy of the semi-detached scenario is restricted by the 
value of the radiative efficiency of accretion K which defines the 
quantity Lacc^K exerting the radiation pressure of a hot spot. 

RLOF 

- — that is available after re- 



Lacc is that part of L"^^ - Gx 
duction due to the fact that matter impinging on the gainer starts 
at the first Lagrangian point and not at infinity. This accretion lu- 
minosity is weakened by the fact that only a fraction can be con- 
verted into radiation and strengthened because the energy of the 
impacting material is concentrated in a hot spot which is signif- 
icantly smaller than the entire gainer's surface. Van Rensbergen 
et al. ([2008) defined a quantity K which enables the calculation 
of the contribution of the hot spot to the total luminosity of the 
gainer It is easier to visualize the action of the hot spot using 
as given by relation (|IJi: 



K = 



1 

K 



spot 



T\f ] 



X (5770)4 



(1) 



Unfortunately, there are only 1 1 reliable hot spot tempera- 
tures available in the literature. Eight systems (VW Cep, CN 
And, KZ Pav, V361 Lyr, RT Scl, U Cep, U Sge, and SV Cen) ai-e 
direct impact systems, whereas three of them (SW Cyg, V356 
Sgr and /3 Lyr) have a transient accretion disk. In the case of the 
formation of a hot spot on the edge of an accretion disk we have 
to replace in equation ([T]): Rg by RiUsk and T^ff^g by Tedge,disk- 
Our liberal evolutionary calculations have thus been gauged with 



small number statistics, producing an empirical relation for the 
radiative efficiency of accretion K concentrated in a hot spot, 
increasing with the total mass of the system as: 



M,, Mg 1-^35 
= 4.386 x[—-H—] 



(2) 



Radiative efficiency, rf, is usually defined through the lumi- 
nosity Ladd which is added to a system as a consequence of the 
transfer of matter at a rate M: 



-'add 



T] M c 



(3) 



A numerical value can hence be calculated from relation (|2| 
for the radiative efficiency of mass accretion by a Main Sequence 
gainer: 



1.735 



77 = 1.88 X 10"^ X [— + — ] xDxS 



(4) 



The factor D is the geometric factor taking into account that 
matter does not fall onto the gainer from infinity but from the 
first Lagarangian point. This factor is zero for a contact system 
and goes to unity as Li goes to infinity. The factor S is the frac- 
tional surface area of the hot spot. Consequently we can compare 
the values of i] calculated with relation (|4| with values that are 
well known from fundamental physics: 77 = 0.007 for complete 
nuclear fusion of hydrogen, 77 grows smoothly from 0.057 for 
a mass gaining non-rotating black hole to 0.32 for a black hole 
rotating at maximum plausible spin (Thome ri974l l. 

3.3. Tidal interaction 

All the systems in the grid have been calculated both with strong 
and weak tidal interaction. The formalism for the tidal interac- 
tion was taken from Zahn d 19771 1, who gives a suitable approxi- 
mation for the synchronisation time-scale: 



Tsync (jr) ^ q X [ — ] 



(5) 



This expression uses the semi major axis a of the binary and 
a mass-ratio q, in which the star that has to be synchronized is in 
the denominator This is the gainer in our case, so that q - 

Tidal interactions modulate the angular velocity of the gainer 
ojg with the angular velocity a>orb of the system. According to 
Tassoul (12000b one can write: 



1 



1 



1 



dt 



X fsyn, 



(6) 



Tidal interactions spin the gainer down when (Dg > cjort- 
Tides spin the gainer up when Ug < cjorh- fsync - 1 represents 
weak tidal interactions whereas fsy„c = 0.1 implies strong tides. 

Strong tidal coupling should be preferable when the spheri- 
cal shape of the gainer is severely elongated due to rapid rotation. 
As a result we found that binaries may lose matter during short 
eras of rapid RLOF which occur soon after the onset of RLOF 
during hydrogen core burning of the donor and sometimes also 
short after the subsequent onset of RLOF during hydrogen shell 
burning of the donor. Binaries with initial primary masses below 
6 Mq certainly show spinning up of the gainer and increased ac- 
cretion luminosities concentrated in hot spots, but the combined 
energies of both events are never sufficient to overcome the bind- 
ing energy of the system: those systems evolve conservatively. 
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Following increasing values of the initial primary mass we find 
that binaries with initial primary masses below 8Mq do not lose 
a significant fraction of the transferred mass. Figures[T]to|3]show 
the amounts of mass lost by systems with initial primary masses 
from 8Mq on, respectively for ^= 0.4, 0.6 and 0.9. Up to 9Mq 
are lost by systems with a 15Mq primary mass at birth. There 
are no huge differences in mass loss between similar cases cal- 
culated respectively with weak and strong tidal interaction. 

4. The distribution of mass-ratios and orbital 
periods of Algols 

4.1. The observed distributions 

Van Rensbergen et al. (12008b used the mass-ratio and orbital 
period distribution of Algols out of a sample of 303 observed 
systems. These systems have been taken from the catalogue of 
Budding et al. (2004) extended with semi-detached Algols from 
the catalogue of Brancewicz et al. ( 119961 ). All these systems are 
issued from a binary with a B-type primary at birth. A majority 
of 268 systems have a late B-type primary progenitor, leaving 
only 35 systems with an early B-type primary at birth. The initial 
conditions described in section ill use q = ^ as a definition of 
the mass-ratio, but since the gainer has become the most massive 
component of the Algol-system, we use q = ^ as the definition 

of mass-ratio of an Algol-system. The observed Algols combine 
a large fraction of systems where Algol characteristics are pro- 
duced from hydrogen core burning of the donor on (case A and 
case A/B) and a smaller fraction where Algols are produced only 
after the ignition of hydrogen shell burning in the donor (case B). 

4.2. Tlie calculated distributions 

Our conservative and liberal binary calculations consider a bi- 
nary to be an Algol when the semi-detached system shows the 
typical characteristics as mentioned by Peters (120011) : 

• The less massive donor fills its Roche Lobe. 

• The most massive gainer does not fill its Roche lobe and is 
still on the Main Sequence 

• The donor is the cooler, the fainter and the larger star 

The calculations in this paper do not include the Case B- 
calculations, so that the final comparison between observation 
and theory will be delayed until the Case B-calculations are also 
included in the CDS-catalogue. The calculations in this paper 
however compare the results as obtained with the assumption 
of strong tidal interaction with those obtained with weak tides. 
Although significant differences frequently occur between both 
assumptions, they yield a very similar overall distribution of 
mass-ratios and orbital periods of Algols. 

4.3. The mass-ratio distribution of Algols with a late B-type 
primary at birth 

Since binaries with a late B-type primary at birth hardly lose 
any mass during their evolution, the mass-ratio distribution as 
obtained with the liberal assumptions will not differ very much 
from the one obtained using the conservative binary evolution- 
ary code. The fact that » 45 % out of a sample of 268 Algols are 
observed with high mass-ratios (q e [0.4-1]) is thus neither re- 
produced with the liberal binary evolutionary scenario nor with 
the conservative code. All these codes produce only » 12 % of 



Algols with a mass-ratio q above 0.4. Figure |4] shows the mass- 
ratio distribution of Algols with a late B-type primary at birth as 
obtained from theory with the conservative and the two different 
liberal assumptions about the strength of the tidal interaction. 

4.4. The mass-ratio distribution of Algols with an early B-type 
primary at birth 

In section (|3]l we showed that binaries with an early B-type pri- 
mary at birth lose a large fraction of the transferred mass dur- 
ing their evolution. The short liberal era during which mass is 
lost occurs when the mass transfer rate is large, i.e. soon after 
the onset of RLOF when the binary is in its pre- or early-Algol 
stage. The mass-ratio distribution as obtained with the liberal 
assumptions will now differ very much from the one obtained 
using the conservative binary evolutionary code. The fact that 
^ 46 % out of a sample of (only) 35 Algols are observed with 
high mass-ratios {q e [0.4-1]) is reproduced with the liberal bi- 
nary evolutionary scenario only 39 %). The large fraction of 
Algols observed with very high mass-ratios 17 % with q e 
[0.85-1]) remains however very hard to explain. Since the liberal 
era occurs mainly in the very early Algol-stage (when q x 1) it 
is clear that the Algol will populate the lower values of q dur- 
ing the much longer lasting eras of quiet and slow mass transfer. 
The binary loses mass during the short lasting era of rapid mass 
transfer and then shows Algol characteristics during long lasting 
eras of quiet and slow mass transfer. Figure [5] shows the mass- 
ratio distribution of Algols with an early B-type primary at birth 
as obtained from theory with the conservative and the two differ- 
ent liberal assumptions. These conclusions are weakened by the 
fact that they are based on a comparison with only 35 Algols. 

4.5. Period distribution of Algols with a B-type primary at 
birth 

Figure |6] shows the orbital period distribution of Algols with an 
early B-type primary at birth as obtained from theory with the 
conservative and the two different liberal assumptions. We find 
no large differences between the conservative and liberal results. 
It will be very hard to reproduce by theory the Algols observed 
with orbital periods below 1 day (» 9 %) since those systems 
always merge. Algols observed with orbital periods above 15 
days (a: 5 %) will be reproduced by theory when one includes 
those Algols which are only formed after ignition of hydrogen 
in the shell of the donor It was e.g. illustrated for the conser- 
vative case by Van Rensbergen (2003) that these cases B pro- 
duce a large fraction of Algols with large orbital periods. These 
Algols will however not influence the final global distribution 
of mass-ratios and orbital periods of Algols very much, since 
their contribution to the Algol-population will not be very large 
due to the fact that the hydrogen shell burning duration will al- 
ways last much shorter than its previously lived hydrogen core 
burning time. Anticipating on finishing our catalogue of binary 
evolutionary calculations with the cases B of RLOF, we may 
conclude that theory reproduces the observations well. 



5. Conclusions 

Mass impinging from the donor spins the gainer up and creates a 
hot spot in its equatorial zone or on the edge of its accretion disk. 
The combined energy of the enhanced rotation and increased ra- 
diation from the hot spot may exceed the binding energy of the 
system. We present the results of liberal evolutionary calcula- 
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Fig. 1. Amount of mass lost by a binary with a B-type primary at birth and an initial mass-ratio jj-^ - 0.4. The amount of mass (in 
Mq) is on the vertical axis. The initial period is on the horizontal axis. The initial mass of the donor and the type of tidal interaction 
are mentioned in the legend. Donors with initial masses below 8 lose only a little amount of mass and are not included. 





tions for binaries with a B-type primary at birth and with small 
initial orbital periods so that RLOF starts during hydrogen core 
burning of the donor. A significant fraction of the transferred 
mass is lost by the system in the case of binaries with an early 
B-type primary at birth. Systems with a late B-type primary at 
birth hardly lose any matter. 



The observed distribution of orbital periods of Algols is well 
reproduced by conservative as well by liberal theoretical calcu- 
lations. 

Although the theoretically calculated liberal mass-ratio dis- 
tribution of Algols with an initial early B-type primary fits the 
observations much better, the overall observed mass-ratio distri- 
bution of Algols still shows too many systems with large mass- 
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Fig. 4. Simulated distribution of mass-ratios of Algols issued from a binary with a late B-type primary at birth and an initial 
orbital period so that Case A RLOF occurs. Since these binaries do not lose a large amount of mass, the mass-ratio distribution is 
approximately the same for conservative evolution as compared to liberal evolution, producing only a; 12 % of Algols with high 
mass-ratios: q e [0.4-1]. 
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Fig. 5. Simulated distribution of mass-ratios of Algols issued from a binary with an early B-type primary at birth and initial orbital 
periods so that Case A RLOF occurs. These binaries lose a large amount of mass so that the mass-ratio distribution differs for 
conservative evolution as compared to liberal evolution, producing ^ 39 % of Algols with high mass-ratios {q € [0.4-1]) in the case 
of hberal evolution only. 
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Fig. 6. Simulated distribution of orbital periods of Algols issued from a binary with a B-type primary at birth and initial orbital 
periods so that Case A RLOF occurs. Conservative and liberal evolution produce almost the same result. 



ratios. This is due to the fact that the calculated era of rapid mass with large mass-ratios. Binary evolutionary calculations yielding 
transfer is very short so that the binary rushes through the states eras of rapid mass transfer lasting for a longer time with some- 
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what lower peak values of the mass transfer rate would produce 
more Algols with large mass-ratios but have never been pub- 
lished. We have compared the mass transfer rates as obtained 
by our binary evolutionary code in the conservative mode with 
mass transfer rates for conservative evolution as produced by 
previous authors. According to Kippenhahn et al. (I1967al lb]l a 9 
Mq donor transfers more than 5 M© to his initial 5 M© gainer 
in 6 X lO'* years during hydrogen core burning of the donor 
and almost 7 Mq in 4 x 10"* years when RLOF starts after ex- 
haustion of hydrogen in the core of the donor A 2 Mq donor 
transfers 0.45 Mq to his 1 Mq companion in 3.1 x 10^ years in 
the rapid phase of mass transfer during hydrogen core burning 
of the donor Paczyhski et al. (1967a b) calculated the conser- 
vative evolution of a binary with a 16 Mq future donor at birth 
and a 10.67 Mq companion. With an initial orbital period lead- 
ing to case A RLOF he finds that almost 8 Mq are transferred 
to the gainer in 4 x lO'* years. When RLOF starts after exhaus- 
tion of hydrogen in the core of the donor a short era of mass of 
mass transfer is found with a peak value as high as 3.4 x 10 

Our calculated durations of rapid mass transfer are very 
similar to those mentioned above whereas our peak values are 
somewhat lower. This is due to the fact that our stellar models 
are calculated with Rogers- Iglesias opacities (T99T) which were 
not available previously. Therefore, our calculated durations and 
peak values of mass transfer rates agree very well with those as 
published by Nelson and Eggleton (120011) for a representative set 
of interacting binaries. The occurrence of many observed Algols 
with large mass-ratios thus remains unexplained. Future inves- 
tigations should explore other interactions between the gravita- 
tional RLOF and the internal thermal structure driving the evo- 
lution of the radius of the donor. 

Podsiadlowki et al. ( |1992| l pointed out that liberal theoretical 
calculations depend very much on the amount of mass lost from 
the system (characterized by the parameter /3) and the amount 
of angular momentum taken away by this matter (character- 
ized by the parameter a). Our liberal code calculates /3{t) self- 
consistently within the model and assumes that matter is lost 
from the hot spot on the gainer (or edge of its accretion disk) so 
that the escaping matter takes only the angular momentum of the 
gainer's orbit. It is clear that if matter would escape at another lo- 
cation (another choice of the parameter a, e.g. characteristic for 
L2 as the position of mass loss from the system) the calculated 
population of Algols could be different. 
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